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PARAMETRIC  STUDY  OF  THE  DYNAMIC  STABILITY  OF  TOWED  VESSELS 

Several  accidents  in  towing  operations  of  barges  or 
disabled  ships  in  restricted  and  open  waters  have  made 
necessary  the  investigation  of  the  course  keeping  stability 
of  towed  vessels.  In  this  work  a  non-linear  mathematical 
model  is  used  to  simulate  the  slow  surge,  sway,  and  yaw 
motions  of  a  vessel  towed  by  a  heavy  catenary  towline.  The 
effect  of  geometric  parameters  of  the  system  on  the 
stability  of  equilibrium  configurations  is  analyzed. 

It  is  shown  that  for  certain  choices  of  towing  system 
parameters,  dynamic  loss  of  stability  may  occur  which 
results  in  qualitatively  different  asymptotic  response.  The 
results  of  this  study  identify  regions  in  the  parameter 
space  that  lead  to  either  safe  operations  or  hazardous 

j  ,  >  -v  ' J'  o. •  •  "T  /'  •'  '-I  •  'j: 

system  response. 
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I.  INTRODUCTION 


A.  BACKGROUND 

A  long  history  of  towing  accidents  resulting  in  loss  of 
life,  damage  to  property,  and  pollution  of  the  environment 
have  prompted  many  studies  into  the  dynamics  of  towing 
operations.  Of  primary  concern  was  the  motions  of  the  towed 
vessel  in  the  horizontal  plane  (yaw,  sway,  and  surge) . 
Excessive  and  unstable  motions  could  lead  to  collisions  and 
capsizing.  The  ability  to  predict  the  motion  of  a 
particular  towing  system  would  be  of  particular  benefit  to 
ship  designers  and  towing  operators,  by  identifying  those 
situations  where  the  towing  operations  would  be  the  safest, 
or  those  which  must  be  avoided. 

Previous  studies  at  the  University  of  Michigan  and 
elsewhere  had  developed  mathematical  models  and  numerical 
techniques  to  analyze  towing  dynamics,  and  had  identified 
those  parameters  which  are  of  primary  importance  to  the 
stability  of  the  towing  system.  The  linear  model  usually 
used  to  describe  ship  motions  [Ref.  1,  Chapter  7]  is 
inadequate  for  the  towing  problem.  Non-linear  models,  as  in 
[Ref.  2],  must  be  utilized  to  accurately  describe  the 
towing  system.  These  studies  had  identified  the  position  of 
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the  towline  attachment  point  on  the  towed  vessel  and  the 

towline  tension  as  the  most  significant  (and  controllable) 

parameters  of  the  towing  system.  • 

B.  PROBLEM  CONDITIONS 

In  this  study,  computer  programs  developed  in  [Ref.  3] 
were  used  to  analyze  the  effect  of  different  parameter 
combinations  on  the  towed  stability  of  three  vessels.  These 
programs  use  bifurcation  to  identify  the  unstable  and 
stable  regions  of  the  parameter  space.  The  principal 
parameters  studied  were  (Fig.l): 

1.  longitudinal  position  of  the  towline  attachment  point 
forward  of  the  towed  vessel's  center  of  gravity,  xp; 

2.  athwartships  position  of  the  towline  attachment  point 
port  or  starboard  of  the  towed  vessel's  centerline, 

yp; 

* 

3.  length  of  the  towline,  Lw. 

In  the  model  used  in  this  study,  unlike  [Ref.  2],  the  * 

towline  is  modeled  as  an  inextensible  caternary,  thus 
making  towline  tension  a  function  of  its  length.  The  model 
conditions  were: 

1.  speed  of  towing  vessel  of  2  knots; 

2.  towing  vessel  on  steady  course; 

3.  calm  seas,  no  wind;  i.e.,  no  external  environmental 
forces . 

Characteristics  of  the  towed  vessel  were  inputted  into  the 
programs  from  a  data  file  containing  hydrodynamic 
coefficients,  resistance  data,  towline  characteristics,  and 
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skeg  data,  if  applicable.  The  effect  of  asymmetrical  forces 
acting  on  the  towed  vessel,  such  as  the  presence  of  a 
propel lor  or  an  environmental  force,  are  introduced  through 
a  bias  in  the  data  file.  All  dimensions  are 
nondimensionalized  with  respect  to  the  towed  ship's  length 
between  perpendiculars  (LBP) . 

Three  vessels  were  studied  (Fig. 2): 

1.  a  191  foot  barge  with  a  skeg,  with  no  propellor 
(i.e. ,  no  bias) ; 

2.  a  1066  foot  tanker  with  no  skeg,  but  with  a  propellor 
(i.e.,  with  a  bias); 

3.  the  same  barge  as  in  1),  but  without  the  skeg  and 
with  a  propellor  (i.e.,  with  a  bias). 

Unlike  previous  studies,  this  work  includes  the  effect 
of  athwartship  position  of  the  towline  attachment  point  in 
the  stability  of  the  towing  system.  * 

Chapter  II  provides  background  into  the  problem 
formulation  and  stability  analysis  used  in  this  study. 

Chapter  III  presents  the  results  of  the  analysis  and 
discusses  some  practical  aspects  of  these  results.  Chapter 
IV  discusses  the  conclusions  which  can  be  made  from  the 
results  of  this  study  on  the  stability  of  the  towing  system 
and  the  use  of  the  techniques  used  herein. 
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II.  PROBLEM  FORMULATION  AND  METHOD  OF  APPROACH 

Slow  motions  of  a  towed  vessel  in  the  horizontal  plane 
are  described  by  a  system  of  six  nonlinear,  coupled, 
differential  equations.  [Ref.  4  and  5]  In  its  standard  form 
this  system  is 
1 

X1  =  ~"tFl(xl,x2'x3)  +  Tsurge(x4'x5'x6) 1 

m  -  *u 

Jz  "  Nr 

x2  _  — - — — — ' [F2 (xi,x2 »x3)  +  Tsway (x4 ' x5 ' x6^ 

D 

Yr 

+  [ f3 (X1 ' x2 ' x3 )  +  xpTsway  ( x4 ' x5 ' x6 )  ^p^surge (x4 ' x5 ' x6^ ' 


x3 


Nv 

[ F2 (X1 ' x2 ' x3 )  +  xpTsway ^X4'X5,X6^  3 


m  -  Y. 


[F3(x 


Xo ,X3) +X 


1'  *2 


pTsway (x4 ' 


x5  /  x6 


) -ypTsurge  ^x4 


(x5,x6) ] 


* 


* 


x4  =  xlCosx6  -  x2sinx6  -  U, 


x5  =  xls^-nx6  +  x2cosx6' 


x6  =  x3' 
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where 


» 
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Tsurge  ( x4 ' x5 ' x6 )  ~  ^x(x4,x5,x6)cosx6+^y(x4'x5' 

_Tsway<x4'x5'x6)  =  Rx(x4'x5'x6>sinx6_Ry(x4'x5' 
D  denotes  the  known  quantity 


x6)sinx6, 

x6>cosx6' 


D  =  (m-Yv) (I2“Nr) -YrNy, 


and 


R1  (xl  >  x2  '  x3 )  ^uxl+^uuxl+^/®^uuuxl+^wx2+^^wux2xl 


2  2 

+  ^xrrx3+?sXrrux3xl+  (xvr+ltl)  x2x3+xrvuxlx2x3  - 


R2  (X1 '  x2  '  x3 )  Y0+Y0uxl+Youuxl+Yvx2+1/ 6Yvwx2+'lYvrrx2x3 
+Yvuxlx2+^Yvuux2xl^  ( Yx'“roX2 )  x3+1/6Yrrrx3  ^^Yrwx3x2 
+Yrux3xi+^ Yruux3xl ' 


2  3  2 

R3  (xi *  x2 ' x3  )  =R0+R0uxl+Rouuxl+Rvx2+^/ 6Nwvx2+^Nvrrx2x3 

2  3  2 

+^vuxlx2+'5^vuux2xl+^rx3't’^-/ ®^rrrx3+'i^rwx3x2 

+Nrux3xl+!*Nruux3xl- 

In  the  above  equations,  x-^  denotes  the  sway  velocity  in 
surge  (longitudinal  motion)  of  the  towed  vessel,  x2  the 
velocity  in  sway  (lateral  motion) ,  x3  the  angular  velocity 
in  yaw  (turning  motion  about  the  vertical  axis) ,  x4  and  x5 
the  coordinates  of  the  center  of  gravity  of  the  towed  vessel 
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with  respect  to  an  (x,y) -coordinate  system  moving  with  the 
towing  vessel,  and  x6  the  towed  vessel  yaw  angle.  Further,  U 
is  the  steady  towing  vessel  velocity  in  the  x-direction,  xp 
and  yp  are  the  coordinates  of  the  towline  connection  point 
on  the  towed  vessel  with  respect  to  an  (X,Y) -coordinate 
system  with  its  origin  at  the  towed  vessel  center  of 
gravity,  and  Rx,  Ry  are  the  towline  restoring  forces.  The 
towing  system  configuration  and  notation  conventions  are 
shown  in  Figure  1.  Expressions  for  Flf  F2,  F3  are  derived  by 
Taylor  expansion  in  terms  of  the  relative  velocities  x-^,  x2, 
x3  of  the  towed  vessel  with  respect  to  the  water.  In 
nonlinear  analysis  terms  up  to  third  order  are  used  whereas 
terms  beyond  third  order  and  second-  and  higher-order 
acceleration  terms  are  usually  neglected.  Subscripts  u,  v,  r 
indicate  derivative  of  force-moment  component  with  respect 
to  xlf  x2,  x3  respectively,  and  subscript  c  indicates 
propellor  dependent  terms.,  which  represent  a  source  of 
system  asymmetry.  These  terms  are  zero  in  the  absence  of  a 
propellor.  Terms  Xabc,  Yabc,  Nabc,  where  a,  b,  c  are  dummy 
independent  variables  representing  u,  v,  r,  are  usually 
called  slow  motion  derivatives.  In  unsteady  reference 
motion,  slow-motion  derivatives  are  considered  as  functions 
of  the  frequency  of  motion.  In  our  study  of  slowly  varying 
reference  motions,  we  assume  that  slow  motion  derivatives 
are  time  independent.  This  is  a  good  approximation  for  ships 
with  usual  hull  shapes  and  moderate  speeds. 
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Rx  and  Ry  denote  restoring  forces  from  the  towline,  and 
for  a  quasistatic  towline  response  they  are  expressed  as 
implicit  functions  of  x4,  x5,  x6.  In  this  study,  the  model 
used  for  the  towline  is  that  of  an  inextensible  heavy 
caternary  with  nonlinear  force-displacement  characteristics 
as  given  in  [Ref.  3]. 

In  compact  notation  the  above  system  of  six  ordinary 
differential  equations  is  denoted  as 

x=f(x)  (1) 

where  x  and  f  are  six  dimensional  vectors.  To  analyze  the 
stability  properties  of  (1) ,  the  first  step  is  to  identify 
the  equilibrium  configuration  of  the  system.  For  this  we 
have  to  solve  a  system  of  six  nonlinear,  coupled  algebraic 
equations 

f (x) =0  (2) 

where  x  denotes  an  equilibrium  configuration.  It  can  be 
shown  [Ref.  5]  that  system  (2)  has  at  most  three  solutions 
in  x  corresponding  to  three  distinct  equilibrium  positions. 
In  this  study  we  concentrated  our  efforts  on  one  of  these 
equilibrium  positions,  namely  the  one  which,  in  the  absence 
of  a  bias  in  the  system,  corresponds  to  the  towed  vessel 
being  located  directly  astern  of  the  tow-tug.  This  is  the 
most  interesting  in  applications.  Having  computed  x,  its 
stability  properties  can  be  established  as  follows: 
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Linearization  of  (1)  around  x  leads  to  the  linear  system 

z=Az  (3) 

where  z  represents  the  excursion  from  the  equilibrium  x,  and  # 

A  is  a  constant  6x6  matrix.  If  all  eigenvalues  of  A  have 

« 

negative  real  parts,  then  x  is  stable;  otherwise  it  is 
unstable. 

In  this  study,  we  performed  parametric  analysis  of  the 
central  equilibrium  in  therm  of  towline  length  Lw,  and  the 
towing  point  coordinates  with  respect  to  the  center  of 
gravity  of  the  towed  vessel,  xp  and  yp.  These  parameters  can 
be  easily  changed  before  or  during  towing  operations  and  can 
provide  a  means  of  passive  control  of  the  towing  system. 

Parameter  Lw  directly  affects  the  amount  of  tension 
developed  by  the  towine.  Parameter  xp  is  directly  related  to 

% 

the  towline  restoring  force  and  moment.  A  small  value  for  xp 

may  not  be  able  to  provide  adequate  restoring  moment  and  may  t 

not  guarentee  system  stability.  On  the  other  hand  a  very 

large  value  for  xp  may  result  in  over-compensation  and 

therefore  instability.  Nonzero  yp  values  result  in  a  source 

of  asymmetry  introduced  in  the  system.  For  a  biased  system 

(for  example  due  to  the  presence  of  a  propellor  or 

environmental  forces) ,  it  should  be  expected  that  an  extra 

appropriate  bias  introduced  via  a  nonzero  yp  helps 

counteract  the  effect  of  the  former  bias,  and  hence  improve 

stability. 
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The  particular  equilibrium  position  will  lose  its 
stability  when  an  eigenvalue  of  the  A  matrix  in  (3)  changes 
its  sign  from  negative  real  part  to  positive  real  part.  The 
case  when  a  real  eigenvalue  crosses  zero  has  ben  analyzed  in 
detail  in  [Ref-  5],  This  corresponds  to  a  static  loss  of 
stability  with  generation  of  additional  equilibrium 
positions  in  the  form  of  solution  branching.  In  this  study 
we  focussed  our  attention  on  the  case  when  a  complex 
conjugate  pair  crosses  the  imaginary  axis.  This  corresponds 
to  a  Hopf  bifurcation:  the  particular  equilibrium 
experiences  a  dynamic  loss  of  stability  and  the  system 
begins  to  oscillate.  The  resulting  periodic  solutions  can  be 
stable  or  unstable,  but  at  any  rate,  such  a  situation  is 
hazardous  and  should  be  avoided  during  towing  operations. 
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111.  RESULTS  AND  DISCUSSION 

A.  BARGE  WITH  SKEG 

The  first  vessel  to  be  studied  was  an  unpowered  barge 
with  a  skeg  aft.  Since  there  is  no  propellor  to  introduce 
a  bias,  the  barge  has  athwartship  symmetry. 

1.  Figure  3:  Critical  Real  Part  vs.  xp 

Program  T0WBIF1  calculates  eigenvalues  for  specific 
Lw  and  yp,  creates  a  file  for  each  of  six  real  and  six 
imaginary  parts,  and  creates  a  separate  file  containing  the 
largest  real  part.  The  real  part  with  the  largest  value  is 
the  critical  indicator  of  the  system's  stability:  if  it  is 
greater  than  zero,  the  system  will  be  unstable;  if  less 
than  zero,  the  system  will  be  stable. 

Figure  3  shows  plots  for  the  critical  parts  for 
yp=0.05  and  three  values  for  Lw.  The  region  where  the  plot 
is  greater  than  zero  indicate  that  range  of  xp  where  the 
system  is  unstable.  For  example,  for  Lw=0.5  the  critical 
real  part  is  greater  than  zero  for  the  range  of  xp=0.18  to 
xp=0 . 48 ,  so  the  system  is  unstable  within  this  range. 

Note  that  as  Lw  increases,  the  unstable  range 
becomes  smaller,  until  for  Lw=3.0  there  is  no  region 
greater  than  zero.  Therefore,  the  system  will  be  stable 
for  all  values  of  xp;  i.e,  the  barge  should  exhibit  no 
unstable  motion. 
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Figure  3.  Critical  Real  Part  vs.  xp  -  Barge  with  skeg 


2.  Figure  4:  yp  vs.  xp,  Lv  as  Parameter 

Program  T0WBIF2  does  the  same  calculations  as 
T0WBIF1  over  a  range  of  values  of  yp  with  a  given  Lw, 
instead  of  a  single  value  of  yp  and  Lw.  In  essence,  Figure 
3  represents  a  cut  of  Figure  4  at  a  single  value  of  yp  and 
Lw.  Unlike  T0WBIF1,  T0WBIF2  writes  a  point  only  where  the 
critical  real  part  changes  sign.  When  plotted,  these 
produce  a  curve  delineating  stable  and  unstable  regions  of 
the  parameter  space.  This  is  the  point  of  the  process;  we 
are  more  interested  in  finding  what  parameters  produce 
stable  or  unstable  system  than  the  actual  results  of  the 
equations  of  motion. 

Recalling  Figure  3,  The  area  inside  the  curves 
represent  the  unstable  region.  For  example,  for  Lw=0.5, 
the  system  is  unstable  for  all  values  of  yp  within  the 
range  xp=0.2  to  xp=0.5.  Increasing  Lw  first  decreases  the 
unstable  range  of  xp  for  high  yp,  then  decreases  the 
unstable  range  of  yp.  For  large  Lw  (>4.0),  the  unstable 
range  virtually  disappears. 

3.  Figure  5:  Lw  vs.  xp,  yp  as  parameter 

Program  T0WBIF3  performs  the  same  calculations  as 
T0WBIF2,  but  with  Lw  as  the  ordinate  and  yp  as  the 
parameter.  Thus  Figure  5  provides  the  same  information  as 
Figure  4  but  with  a  different  perspective. 
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Figure  5.  Lw  vs.  xp,  yp  as  parameter 


As  in  Figure  4,  the  unstable  region  is  inside  the 
curves.  It  clearly  shows  how  increasing  Lw  decreases 

»  the  unstable  range  for  a  constant  yp,  as  was  evident  in 

Figure  4.  It  also  shows  that,  for  constant  Lw  greater  than 
about  0.7,  increasing  yp  also  decreases  the  unstable 
region.  In  the  narrow  range  of  Lw  from  0.2  to  0.7, 
increasing  yp  increases  the  extent  of  the  unstable  range  of 
xp.  This  effect  is  apparent  in  Figure  4,  but  more 
dramatically  presented  in  Figure  5.  It  would  appear  that 
using  two  views  of  the  data  would  emphasize  aspects  of  the 
curve  that  may  be  overlooked  with  one  view. 

Since  positive  values  of  yp  represent  port  side 
placement  of  the  towline  attachment  point,  and  negative 
values  starboard  side  placement,  both  positive  and  negative 

¥ 

values  for  yp  were  studied.  As  expected  from  the  port- 

•  starboard  symmetry  of  the  barge,  curves  for  positive  and 
negative  values  of  yp  were  identical,  and  only  positive 
values  were  presented  here. 

From  an  operational  point  of  view,  one  may  conclude 
from  these  curves  that  for  the  unpowered  barge,  placing  the 
towline  on  an  attachment  point  to  either  side,  as  far 
forward  as  possible,  will  make  the  tow  stable  for  the 
greatest  range  of  towline  length,  but  the  towline  should  be 

*  kept  no  shorter  than  the  length  of  the  barge. 
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B.  TANKER 


The  second  vessel  studied  was  a  tanker  typical  of  those 
now  in  service.  The  effect  of  the  tanker's  propellor  makes 
the  hull  asymmetrical;  this  effect  is  represented  by  a  bias 
included  in  the  tanker  data  file. 

1.  Figure  6:  Critical  Real  part  vs  xp 

Figure  6  plots  data  generated  from  T0WBIF1  with 
yp=0.10  and  two  values  for  Lw.  These  plots  show  the  stable 
region  to  be  between  the  two  zero  values  for  the  curve. 

Note  that  the  stable  region  becomes  smaller  with  increasing 
Lw. 

Note  also  that  both  curves  are  discontinuous  in  their 
slopes.  The  critical  real  parts  file  is  a  composite  of 
several  results  files,  each  of  which  is  critical  over  a 
certain  range.  Each  results  file  forms  a  smooth  curve; 
thus  the  curves  plotted  on  each  T0WBIF1  figure  may  be 
combinations  of  the  critical  section  of  several  results 
files. 

Finally,  note  that  the  stable  region  occurs  over  a 
narrow  range  of  xp,  unlike  the  barge  with  skeg  discussed 
earlier. 

2.  Figure  7:  Lw  vs.  xp,  yp  as  Parameter 

Figure  7  was  produced  from  data  generated  by  T0WBIF3 
for  positive  values  of  yp.  As  was  shown  in  Figure  6,  the 
stable  region  is  inside  the  curves.  The  vertical  line  at 
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Figure  6.  Critical  Real  Part  vs.  xp  -  Tanker 


TANKER  W/  CATENARY 


Figure  7.  Lw  vs.  xp,  yp>0  as  parameter 


xp=0 .404  is  a  common  crossing  point  for  all  curves.  Each 
curve  is  formed  by  two  cusps,  with  the  upper  cusp 
»  dominating  with  decreasing  yp.  Each  cusp  is  the  plot  of 

different  critical  pair  of  eigenvalues;  the  "nose"  in  the 
curves  is  the  point  where  they  intersect. 

Note  how  the  stable  region  gets  smaller  with 
decreasing  yp  for  Lw  less  than  0.7,  for  example,  with 
yp=0 . 0  and  Lw  greater  than  1.0,  there  is  a  very  narrow 
range  of  xp  where  stability  can  be  assured. 

3.  Figure  8:  Critical  Real  Part  vs.  xp 

TOWBIF1  was  again  used  to  form  Figure  8,  this  time 
with  one  value  for  Lw  (Lw=0.6)  and  three  negative  values 
for  yp.  The  negative  yp  curves  pass  frcm  stanle  to  unstable 
regions,  with  the  stable  ranges  for  xp  getting  smaller  as 

w 

xp  becomes  more  negative. 

As  in  Figure  6,  the  curves  are  composites  of  those 
results  curves  which  are  critical  over  a  particular  range 
of  xp. 

4.  Figure  9:  Lw  vs.  xp ,  Negative  Values  of  yp  as 
Parameters. 

Figure  9  data  was  generated  from  TOWBIF3,  with  yp=0.0 
curve  included  to  provide  continuity  with  Figure  7 . 

The  stable  region  gets  smaller  as  yp  decreases  from 
'  0.0.  At  yp=-0.10  the  "nose"  between  upper  and  lower  cusps 

appears  to  be  tipping  up,  with  the  region  inside  the  "nose” 
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Figure  9.  Lw  vs.  xp,  yp<0  as  parameter 


being  stable.  The  curves  plotted  in  Figure  8  were  formed 
using  a  value  of  Lw  which  cut  through  this  nose,  thus 

forming  the  sinuous  curves  which  pass  in  and  out  of  the  < 

stable  region.  Note  that  for  the  most  negative  values  of 
yp,  the  cusps  have  disappeared,  and  the  stable  range  of  xp 
is  slightly  increasing. 

5.  Figure  10:  Lv  vs  xp,  Negative  Values  of  yp  as 
Parameters . 

Figure  10  is  a  "close-up"  of  Figure  9,  focusing  on 
what  is  happening  around  yp=-0.10.  The  lower  cusp  tips  up 
and  merges  with  the  upper  to  form  a  single  curve.  Note  how 
the  stable  range  of  xp  virtually  disappears  for  Lw  greater 
than  0.5  for  yp=-0.10  and  -0.11.  As  was  seen  in  Figure  9, 
stable  range  for  xp  for  Lw  greater  than  0.5  reappears  with 

* 

yp  less  than  -0.15. 

4 

C.  BARGE  WITHOUT  8 KEG 

The  third  vessel  was  a  self-propelled  version  of  the 
barge  studied  in  Section  A  (not  under  power  during  tow) , 
but  without  the  skeg.  As  with  the  tanker,  the  presence  of 
the  propellor,  simulated  by  a  bias  in  the  data  file, 
introduces  port-starboard  asymmetry. 

l.  Figure  11:  Critical  Real  Part  vs  xp  for  Lv=1.5 
Figure  11  shows  curves  for  three  values  of  yp 
(greater  than  zero,  zero,  and  less  than  zero)  and  one  value 
of  Lw.  The  curves  show  the  stable  range  of  xp  to  be  between 
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Figure  10.  Lw  vs.  xp,  yp<0,  Close-up 


the  zero  crossing  points  of  the  critical  real  parts,  as  in 
the  tanker  case.  Also  similarly  to  the  tanker,  the  stable 
region  increases  with  increasing  yp.  These  results  are 
opposite  to  the  propellor-less  barge  with  the  skeg. 

2.  Figure  12:  Lw  vs  zp 

Figure  12  dramatically  shows  how  decreasing  yp 
reduces  the  stable  region.  The  vertical  line  at  xp=0.14  was 
common  to  all  values  of  yp  greater  than  and  equal  to  zero. 
For  values  of  yp  less  than  zero,  the  smooth  shape  of  the 
curve  is  apparent. 

The  tanker  and  self-propelled  barge  cases 
dramatically  demonstrate  the  effect  that  a  bias,  like  a 
propellor,  can  introduce  to  the  stability  of  the  system. 

D.  PRACTICAL  OBSERVATIONS 

Analysis  of  the  graphs  suggests  some  general  principles 
which  may  be  applied  when  conducting  slow  speed  towing 
operations  with  the  vessels  discussed  in  this  chapter. 

While  these  principles  are  of  course  not  generally 
applicable  to  all  vessels,  they  illustrate  how  the  analysis 
techniques  employed  in  this  work  can  be  applied  to  other 
vessels. 

For  the  unpowered,  symmetric  barge  with  a  skeg,  the 
operator  should  have  the  towline  attachment  point  as  far 
out  to  either  side  as  possible  and  the  towline  as  long  as 
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Figure  12.  Lw  vs.  xp,  yp  as  parameter 


practical.  The  attachment  point  can  then  be  placed  at  any 
location  forward  of  the  center  of  gravity  with  stability 
)  assured.  Conversely,  if  the  attachment  point  must  be  on 

the  centerline,  placing  it  as  far  forward  as  possible 
(about  half  the  barge's  length  forward  of  the  center  of 
gravity)  will  assure  stability  for  all  towline  lengths. 

For  vessels  with  an  asymmetrical  bias  (e.g.,  with  a 
propellor) ,  but  without  skegs,  the  attachment  point  needs 
to  be  as  far  to  the  biased  side  as  possible  (in  the  cases 
of  the  tanker  and  self-propelled  barge,  the  +yp  or  port 
side)  and  placed  forward  of  the  center  of  gravity  the 
distance  indicated  on  the  graph  for  all  towline  lengths. 
Placement  of  the  attachment  point  on  the  opposite  side  (in 
the  cases  studied,  the  starboard  side)  will  virtually 
assure  the  system  to  be  unstable, 
ft 


29 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

This  study  highlighted  the  effect  of  athwartship 
position  of  the  towline  attachment  point.  The  common 
assumption  among  ship  operators  prior  to  this  reseach  held 
that  placing  the  towline  on  the  centerline  on  the  foremost 
point  of  the  towed  vessel  would  create  the  most 
advantageous  towing  situation.  Studies  such  as  [Ref.  2] 
have  shown  that  towing  stability  can  be  dependent  on  the 
longitudinal  placement  of  the  attachment  point.  This 
research  has  shown  that  for  certain  conditions,  attaching 
the  towline  off  the  centerline  can  also  improve  towed 
stability.  The  optimum  towing  configuration  requires  a 
combination  of  all  three  parameters  -  longitudinal  and 
athwartship  placement  of  the  towline  attachment  point,  and 
towline  length. 

The  bifurcation  technique  used  in  this  study  can  be 
used  to  produce  stability  information  useful  to  ship 
designers  and  towing  operators.  Stability  information  can 
be  assembled  into  a  convenient  graphical  form  that  clearly 
defines  the  regions  of  stable  and  unstable  operation  based 
on  the  parameters  the  operator  has  the  most  control  over  - 
the  placement  and  length  of  the  towline. 


30 


For  the  ship  designer,  this  technique  can  be  useful  in 
determining  the  implications  particular  design  decisions 
would  have  on  the  vessels  performance  under  tow. 

Depending  on  the  vessel's  use,  adjustments  to  the  design 
can  be  made  to  improve  towing  stability,  or  the  customer 
can  be  forewarned  to  avoid  certain  kinds  of  operations. 
Since  nearly  all  vessels  are  towed  at  some  time,  towing 
performance  should  be  analyzed  for  all  vessels. 

For  the  towing  operator,  this  technique  can  provide 
readily  available  information  about  how  a  particular  vessel 
will  respond  under  tow.  The  operator  can  then  adjust  the 
towing  parameters  ^.,  placement  the  attachment  point 
and/or  length  ci  ,ne  towline)  so  the  tow  will  be  in  its 
most  stable  condition,  or,  if  unavoidable,  know  that  a 
particular  towing  situation  will  be  potentially  dangerous 
and  make  preparations  to  deal  with  it. 

Since  ship  data  is  inputted  through  a  data  file,  the 
towed  performance  of  any  vessel  can  be  analyzed  with  this 
method,  including  structures  such  as  offshore  oil 
platforms.  Existing  vessels  can  be  analyzed,  as  well  as 
different  loading  conditions. 
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Two  principal  disadvantages  are  associated  with  this 


technique: 

1.  The  programs  are  dependent  on  the  quality  of  the  ( 

data  provided.  Determining  hydrodynamic  coefficients 

and  resistance  data  requires  tow  tank  experiments  and 
analysis,  and  are  not  obtained  for  most  vessels; 

2.  The  programs  require  large  amounts  of  computer  time 
and  memory  to  run,  which  may  not  be  available  or  too 
costly  for  potential  users,  especially  to  run 
extensive  "what-if"  scenarios.  This  problem  may  be 
alleviated  as  more  inexpensive,  high  speed,  high 
capacity  micro-  and  personal  computers  become 
available. 


B.  RECOMMENDATIONS 

This  study  was  done  for  only  one  set  of  conditions. 

Further  research  can  be  done  in  determining  the  effect  of 
varying  conditions,  such  as  different  speeds  or  maneuvering 
by  the  towing  vessel,  on  the  stability  of  the  tow.  External 

% 

forces  are  modelled  by  the  bias  in  the  data  file.  A 

systematized  method  of  introducing  biases  into  the  data  ( 

would  enable  the  analysis  of  the  effect  of  environmental 
conditions  on  the  towing  system. 

Further  work  should  be  conducted  to  improve  the  "user- 
friendliness”  of  the  programs.  As  currently  configured,  the 
programs  must  be  run  instructively,  and  graphics  produced 
offline.  This  is  a  time  consuming  process  which  does  not 
use  the  full  capabilities  of  either  the  programs  or  the 
graphics  capabilities  of  the  mainframe.  Program  * 
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improvements  should  focus  on  streamlining  computations  and 
user  interaction,  and  incorporating  graphics,  with  the  goal 
of  making  it  available  as  a  ship  design  tool. 


I 
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APPENDIX 


Driver  programs  used  in  this  thesis  are  shown  here. 
Subroutines  can  be  obtained  by  contacting: 

Prof.  F.A.  Papoulias 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  CA  93943 

FILE:  T0WBIF1  FORTRAN  A1 


C 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 


c 


PROGRAM  T0WBI1 

BIFURCATION  ANALYSIS  OF  TOWING  SYSTEMS 
PARAMETER  DEPENDS  ON  IPAR 
IP AR  *  I  :  XP 

2  :  YP 

3  :  LW 

IT  NEEDS  SUBROUTINES  FROM  TOWING. FTN 
IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  MAS3P • NVD . NV . NRD . I ZZ , NR . LB .LEN ,LW , 
l  NO.NOU.NOUU.NVVV.NVRR -NVDD.NVU.NVUU.NRRR ,NPVV. 

>  NRDD • NR U  » NR UU . ND • NDDU , NDVV , NDRR . NDU.NDUU . NVRD 

DIMENSION  IV1(6) .A<6.6).VV<5>.Xt6).WR(6).WI (6 ) .2(6 .6 ) .SV2 C6 ) 

CCMMON/INTGR/ISKEG.NREDP. ITYS. ID, IFDS. ISTAB, IPROP 
COMMCN/SPAR/MASSP .LW.XPP.YPP.LB 
CC“MCN/SURGE/SUC7) 

common/xsurg/xu  .  XUU .  XU'JU 

CCM*CN/SWAY/SW( 25 ) 

CCMMCN/YAW/YAI 16 ) 

COMMQN/MTER/VCAR .RHO. AB3 -CONI .CON2 
COMMON/RE3IST/VELC40) .RESI(AO) 

COMMON/VELE/UEL( 100 ) 

C0*^MCN/P03TN/XI  ,  Y I  ,  Z I 
CCMMCN/GEOM/AL .RW.G. AET .hw.hwi 
CCMMCN/PROP/ALE-P.EY.DIA.ANIU 
COMVCN/CTNR/XC ( ?9 ) . YC  C  99 ) . TC( 99 1 . TC ( «• ) 

CCMMON/INT1/IC 
COMMCN/DOC/UC .ALPHA 
C0MMCN/UEPT/RLX.RLY.RL2 

CCMMON/SLOPE/PDRXX.PDRXY.PDRYX.PDRYY 
C0KM0N/SLAN/RXX6 .RYY6.RXX.RYY 

OPEN  (UNIT*3$,FILE*‘ BARGES '  ,STATU$«,OLD' ) 

OPEN  ( UN I T *  1 .FILE*'RES0’ .STATUS* ’ NEW'  ) 

OPEN  (UNIT* 11 .FILE*' RESIR ' .STATUS* *  NEW’  ) 

OPEN  (UNIT*12,FILE* 'RESCR* .STATUS* *  NEW ) 

OPEN  ( UNI T* 1 1 -F IlE* ' RES3R’ .STATUS* 1  NEW'  ) 

OPEN  ( UN I  T  *  1 4 .FILE** RE SAP1 .STATUS* 'NEW ) 

OPEN  (UNIT-15. FILE*‘RES5R’ .STATUS* 1 NEW ) 

OPEN  (UNIT*16.FILE*'RES6R' .STATUS* *  NEW > 


TOMOOOIO 
TOWOOCSO 
TOWOOOJO 
TOWOOOAO 
TOW00050 
TOU00060 
TOW00070 
TOWOOCBO 
TOW00090 
TOM00100 
TOWOOl 10 
TOW00120 
TOWOOl 30 
TOWOOl AO 
TOW00150 
TOWOOl 60 
TCW001 70 
TOWOO10O 
TOWOOl 90 
TOW00200 
TOW00210 
TOW00220 
TOWOOCIO 
TOW00ZA0 
TOW00Z50 
TOW00260 
TOW00270 
TOW00280 
TOW00290 
TOWOOSOO 
TOW00310 
TOW00320 
TOW00330 
TOW003A0 
TOW00350 
TOW00360 
TOW00S70 
TOW00380 
TOW00390 
TOWOOAOO 

towooaio 

TOW00A20 

TOW00430 

TOWOOAAO 
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IP  (IPar.EQ.2)  YPP=aa 
IF  (IPAR.E0.2J  LW  * AA 
Al  r  .W-LB-C .20A8D0 

ale=al 

CALL  STABIL( IVV.VV. ISOl) 

C 

C  SET  V*VV(K)  PCR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 

C 

IF  ( IVV.NE.NEQL  )  GO  TO  I 
V  =  VV( IKE  ) 

I*  (CASOCVJ.GT.  1  .DO)  STOP  11)1 
CALL  ESUlLB'V.X.RES-RX.fiY) 

CALL  LINEARCX.RES.A.RX.RY J 

CALL  RG(6 .  6 .A. MR. M2 . MATZ.Z « IV1 ,SV2  #  IER  1  > 

IF  ( IER] . n£ . 0  )  STOP  2222 
CALL  DEG3T8 ( CEOS . WR  3 
WRITE  (1.10)  AA.DFCS 
DO  11  Jal.t 
JR=10*J 

WRITE  (JR. 10)  AA.WR(J) 

JI *2C*J 

WRITE  t  J  I  .10)  AA.WKJ) 

11  CONTINUE 
1  CONTINUE 
50C  STOP 

10  FORMAT  (ZD:0 .10) 

1001  FORMAT  c  ENTER  1  i  XP  VARIATION'./, 

1  •  2  :  YP  VARIATION' ,/. 

2  *  Z  :  LM  VARIATION- ) 

IOC:  FORMAT  ('  ENTER  PARAMETER  RANGE’) 

1  OS 2  FORMAT  c  ENTER  NUMBER  OF  INCREMENTS’) 

1005  FORMAT  (*  ENTER  EQUILIBRIUM  NUMBER  * ) 

1 0C6  FORMAT  ('  ENTER  ESTIMATED  NO.  OF  EQUILIBRIA') 

2001  FORMAT  (215) 

E'C 


TOW00720 
TCW007A0 
TCW00750 
TCW00760 
TCWOC770 
7CWC0780 
TCW007«0 
TOWOC8CO 
TCW0C810 
TON008:0 
TCWC0820 
TOW008A0 
TOWC0850 
TOW00860 
TOW00870 
TOW0088 0 
TCW0C890 
TOW00900 
TQW00910 
TCW00920 
TOW00920 
TOWOOBAO 
TOW009SO 
TCW00960 
TCW00970 
TCU0CB80 
TCW00990 
TCW0100Q 
TOWOIOIO 
TCW01020 
TOWOIOIO 
TCW010A0 
TOWOIOSO 
TCW01060 
TCWCJ070 
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FILE:  T0WBIF2  FORTRAN  A) 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 


c 

c 


c 


c 


PROGRAM  TOWBI2 
PROGRAM  TOWBIF.FTN 

BIFURCATION  ANALYSIS  OF  TOWING  SYSTEMS 

PARAMETERS  ARE:  Xp  .  Yp 

IT  NEEDS  SUBROUTINES  FROM  TOWING , FTN 

USER  DEPENDENT  SUBROUTINES: 

DEGSTB  *  CURVES  ENCLOSING  REGION  II  OF  FIGURE  13 
(SUBROUTINE  DEGSTB  IS  IN  TOWING. FTN) 

DS1  *  CURVES  ENCLOSING  REGION  V  OF  FIGURE  13 
(SUBROUTINE  DS1  IS  IN  SPM81F) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  MASSP .NVD.NV.NRD. IZZ.NR.LB.LEN.LW, 

1  NO • NCU . NOUU . NVVV . NVRR >  NVDD . NVU , NVUU . NRRR . NR VV . 

2  NRDD >  NRU  t NRUU >  ND » NDDD • NDVV . NDRR . NDU  • NDUU . NVRD 

DIMENSION  IVl(6)»A(o>6)tVV(3)'X(6)>WR(6)'WI(()»Z(6<6)t S V2 ( 6 ) 

CCMM0N/INTGR/I3KEG.NREDP , ITYS, ID, IFDS. ISTAB. IPROP 
COM  10N/SPAR/MASSP ,LW . XPP , YPP . LB 
CCMMON/SURGc/SUC  7 ) 

COMKCN/XSURG/XU . XUU , XUUU 
COMMCN/SWAY/SW(15) 

COMMON/ YAW/YA (16) 

COMMCN/MTER/ VCAR . RHO . ABS . CONI . CON2 
COMMON/RESIST/VELUO  )  ,RESI  (AC  ) 

COMMON/VELE/UEL (1 00 ) 

ccmmcn/postn/xi .yi ,zi 

COMMON/GEOM/AL.RW.G. AET.HW.HW1 
COMMCN/PROP/ALE.P.EY.DIA.AWJU 
CCMMCN/CTNR/XC ( 99 ) . YC ( 99 ) ,ZC( 99 ) , TCC  99 ) 

COMMON/ 1  NT  1 / 1 C 
COMMON/ DCC/UC • ALPHA 
CCMMON/UE»T/RLX . RLY ,RLZ 

commcn/slope/pdrxx.pdrxy.pdryx.pdryy 

CCMmCN/SLAN/RXX6 • RYY6 -RXX  »RYY 

OPEN  (UNIT*35,FILE*’ BSKEG2 * »STATUS* ’ OLD’ ) 

OPEN  ( UNI T* 11. FILE*’ RESIR ’  , STATUS*  *  NEW*  ) 

OPEN  (UNIT*12,FILE*'RES2R' .STATUS*’ NEW ) 

OPEN  ( UNI T* 1 2 ,FILES ' RESIR ' .STATUS* ’ NEW ) 

OPEN  (UNIT* 1 A . F ILE* ‘RE GAR  * ,STATUS*'NEW ) 

OPEN  (UNIT* 15 -FILE* 'RES5R* .STATUS*’NEW ) 

OPEN  (UNIT* 1 6 ,F ILE* ’RESoR ’ ,STATUS«’NEW ) 

CALL  INPUTDC10) 

VCAR  *VCAR*1 .689D0 

AL  *LW*LB"0.30A8D0 

ALE  *AL 

MATZ  *0 

I  FLOW* 1 

EPS  *l.D-5 

ILMAX* 1500 


TOWOOOIO 
TOMOOOZC 
TOWOOOJC 
TOWOOOAO 
TOW00050 
TOWOOO  6  0 
TOW00070 
TOW00080 
TOW00090 
TOMOOIOO 
TOWOOllO 
TOW00120 
TOWOOllO 
TOWOOIAO 
TOW00150 
TOW 00 160 
TOW00170 
TOW 00 180 
TOW0019Q 
TOWOOCOO 
TOW002 10 
TOW00220 
TOW00220 
TOW002AU 
TOW00250 
TOW00260 
TOW0C270 
TOW00280 
TCW00290 
TOWOOIOO 
TOWOOllO 
TOWOOI20 
TOWOOIIO 
TOWOOIAO 
TOWOOI50 
TOWOOIbO 
TOW00270 
TOW00I80 
TOW00I90 
TOWOOAOO 
TOWOOAIO 
TOWOOACO 
TOWOOASO 
TOWOOAAO 
TOW00A50 
TOWOCAoO 
TOW00A70 
TOM00A80 
TOW00A90 
TOW00500 
TOW00510 
TOW00520 
TOW00530 
TOW005A0 
TOW00550 
TOW00560, 


WRITE 

.  1001) 

TOW00570 

READ 

.■) 

A 1  ,  A2 

TOW00580 

WRITE 

.1002) 

TOW00590 

READ 

,■) 

NUM1 

TOWOObOO 

WRITE 

.1001) 

TOWOOMO 

READ 

.*) 

bi  .b: 

TOWOC620 

WRITE 

•1004) 

TCW00b30 

READ 

NUM2 

TOWOObAO 
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WRITE  (-.IOCS) 

READ  (•*,*}  IKB 

WRITE  (*.1006) 

READ  (*.")  IDS 

DO  1  I*1.NUH1 

WRITE  (*,2001)  I.NUMJ 
YPP=A1*(A2-Al  )»tI-n/tNUMl-n 

XPP=B1 

CALL  STAB  I L  C I VV » VV  *  I SOL  ) 

SET  V*VV(K)  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 
VsVVtlKB) 

if  (dabgivj.gt.i.do)  stop  nil 

CALL  ECUILB(V.X.RES.RX.RY) 

CALL  L INEAR (X. RES. A.RX.RY) 

CALL  RGC6.6.A.WR.WI .MATZ.Z. IV1 *SV2  » IER1 ) 

IF  (IER1 .NE.O)  STOP  2222 
IF  (IDS.EC.l)  CALL  DEGSTBtDEOS.WR) 

IF  (IDS. EC. 2)  CALL  D5UDE0S.UR) 

DE0SCC*DEC3 
XPCO  *XPP 
L  *0 
DO  2  J=2.NUM 2 
WRITE  («.»)  J 

XPP=EWB2-B1  )«U-l)/lNUM2-n 
CALL  STABILC IVV • VV  *  I SOL) 

SET  V*VV(K )  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 
V*VV< IKB) 

IF  (DABS(V).GT.l.DO)  STOP  1111 

call  EOUILB(V.X.RES.RX.RY) 

call  LINEAR(X.RES.A.RX.RY) 

call  RGCb.b.A.WR.WI.MATZ.Z.lVl ► SV2  # IER 1 ) 

IF  (IERl.NE.01  STOP  2222 
IF  (IDS.EO.l)  CALL  DEGSTBCDEOS.WR ) 

IF  (ICC. EC. 2)  CALL  DS1(DE0S«WR) 

DEOSNN=DEOS 

XPNNsXPP 

PR=DEOSOC*DEOCNN 

IF  (PR.GT.O.DO)  GO  TO  3 

L*L*1 

IF  (l.GT.6)  STOP  1 00 0 
I L  -  0 

HP  0  =  XPCO 

XPN*XPNN 

DECSO*DECS03 

DECSNs  DECSNN 

XPL  *  XP  G 

XPPsXPN 

DECSL«CEOSO 

DECSR  *CEOSN 

XPP*(XPL*XPR  )/2 . DO 

CALL  STABILUVV.VV.ISOL) 

V* VV (IKB) 

if  (dabscv).gt.i.do)  stop  nil 

CALL  ECUILBIV.X.RES.RX.RY) 

CALL  LINEAR(X.RES.A.RX.RY) 

CALL  RG(6.b.A.WR.WI,MATZ.2«lVl #SV2 •  IER1 ) 

IF  ( IEP1 .NE.O)  STOP  2Z22 
CALL  DEGSTBtDEOS.WR ) 

DEOSM.DEOS 


TOWOObSO 
TOwOObbO 
TOWOOb  70 
TOWOObBO 
TOW00b90 
TOW00700 
TOW007 1 0 
TOW00720 

TOW00730 
TOW00740 
TOW00750 
TOW007b0 
TOW00770 
TOW00780 
TOW00790 
TOMOOBOO 
TCW00810 
TOW00820 
TOW00830 
TOW008A0 
TOU008SQ 
TOW0Q860 
TOW00870 
TOW00880 
TOWC0890 
TOW00900 
TOW0091O 
TOU00920 
TOW009I0 
TOW00940 
TOW00950 
TOW00960 
TOW009/0 
TOW00980 
TOW00990 
TOWQIOOO 
TOWOIOIO 
TOW0I020 
TOW0I030 
TOWOl 040 
TOW0I050 
TOWOlObO 
TOW01070 
TOW01080 
TOW01090 
TOWOIIOO 
Towomo 
TOWOl 120 
TOWOI1SO 
TOWOl 140 
TOWOl 150 
TOWOl 160 
TOWOl 1 70 
TOWOl 180 
TOWOl 1 90 
TCW01200 
TOW01210 
TOW01220 
TOW01230 
TOWOl 240 
TCW012S0 
TOWOl 2b 0 
TOWOl 270 
TOWOl 280 


XPM*XPP 

PRL*DEOCL"DECSM 

PRR*DEOSR*DEOSM 

IF  tPRL.GT.O.DO)  GO  TO  5 

XPO=XPL 

XPN*XPM 

DEOSO=DEOSL 

DEOSNOEOSM 

IL*IL*1 

IF  (IL.GT.ILHAX)  STOP  3100 
DIF* DABS ( XPL-XPM ) 

IF  (DIF.GT.EPS)  CO  TO  6 

XP=XPH 

GO  TO  4 

5  IF  (PRR.GT.0.D0)  STOP  3200 
XPO»XPM 


TOWO 1290 
TOWO 1300 
TOW01310 
TOWO 1320 
TCW01S30 
TOWOl 340 
TOWOl 350 
TOW01360 
TOWOl 370 
TOWOl 380 
TOWOl 390 
TOWOKOO 
TOWO 14 10 
TOWOl 420 
TOW01430 
TOWOl 440 


XPN«XPR 
DEOSO=DEOSM 
DE03N*DE0SR 
IL* IL* 1 

IF  CIL.GT.ILMAX)  STOP  3100 
DIFsDABGtXPM-XPR) 

IF  (D1F.GT.EPS)  GO  TO  6 
XP«XPM 

4  LLL= 1 0  *L 

WRITE  (LLL.10)  XP.YPP 
3  XPOOsXPNN 

DEOCOO*DEOSNN 

2  continue 

1  CONTINUE 
STOP 

10  FORMAT  (2D20.10) 

1001  FORMAT  C  ENTER  RANGE  OF  Yp  VARIATION*) 

1002  FORMAT  (*  ENTER  NUMBER  OF  INCREMENTS  IN  Yp* ) 

1003  FORMAT  (*  ENTER  RANGE  OF  Xp  VARIATION*) 

1004  FORMAT  (*  ENTER  NUMBER  OF  INCREMENTS  IN  Xp') 

1005  FORMAT  (*  ENTER  EQUILIBRIUM  NUMBER*) 

1006  FORMAT  (•  ENTER  DEGREE  OF  STABILITY  CONTROL*  ) 
2001  FORMAT  (215) 

END 

C 

SUBROUTINE  DSl(DEOS.WR) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-2) 

DIMENSION  WR ( 6 ) 

DEG31«-1 .D30 
DO  1  1*1.6 

IF  (WR( I 1.LT.DEOS1 )  GO  TO  1 
DEOS1 *WR ( I ) 

IJ*I 

1  CONTINUE 
DEOC2*- 1 .D30 
DO  2  1*1.6 

IF  (1J.EO.I)  GO  TO  2 

IF  ( WR ( 1 ) . LT . DEOS2 )  GO  TO  2 

DEOS2*WR ( 1 ) 

!JJ*r 

2  CONTINUE 
DECS*- 1 . D30 
CO  !  1*1.6 

IF  (I.EO.lJ.OR.i.EO.IJJ)  GO  TO  3 
IF  (WR(I).GE.DEOS)  DEOS*WR(I) 

3  CONTINUE 
RETURN 
END 


TOM01450 
TOW01460 
TOW01470 
TOWOl 480 
TOW01490 
TOW01500 
TOW01510 
TOW01520 
TOWO1530 
TOWOl 540 
TOW01550 
TOW01S60 
TOWO 1570 
TOWOl 580 
TOWO 1590 
TOWO 1600 
TOWOl 610 
TOWO 1620 
TOW01630 
TOWO l 640 
TOWOl 650 
TOWO 1660 
TOWOl 6 70 
TOWO 1680 
TOWO 1690 
TOWO 1700 
TOW01710 
TOWO l 720 
TOWOl 730 
TOWO 1740 
TOWO 1750 
TOWOl 760 
TOW01770 
TOWOl 780 
TOWOl 790 
TOW01800 
TOWO 18 10 
TOW01820 
TOWO 18 30 
TOWO 1840 
TOWO 1850 
TOWO 18 60 
TOWO 1870 
TOWO 1880 
TOW01890 
TOWOl 900 
TOWOl 910 
TOW01920 
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FILE:  T0WBIF3  FORTRAN  A1 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 


c 


c 


PROGRAM  TOkBIF3.FTN 

BIFURCATION  ANALYSIS  OF  TOWING  SYSTEMS 

PARAMETERS  ARE:  Xp  .  Lw 

IT  NEEDS  SUBROUTINES  FROM  TOWING. FTN 

USER  DEPENDENT  SUBROUTINES: 

DEG3TB  *  CURVES  ENCLOSING  REGION  II  OF  FIGURE  13 
(SUBROUTINE  DEGSTB  IS  IN  TOWING. FTN) 

DS1  *  CURVES  ENCLOSING  REGION  V  OF  FIGURE  13 
(SUBROUTINE  DS1  IS  IN  SPMBIF) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-D 

DOUBLE  PRECISION  MA3SP .NVD.NV.NRD. 122, NR ,LB .LEN.LW. 

1  NO • NOU . NOUU • NVVV • NVRR . NVDD . NV  U . NVUU . NRRR • NR VV . 

2  NRDD • NR  U . NR  UU . ND . NDDD • NOW . NDRR . NDU . NDUU . NVRD 

DIMENSION  IV1(6)<A(6>6)<VV(3),X(6).WR(6)<W1(6). 2(6.6). SV2 ( 6 ) 

CCMMON/INTGR/ISKEG.NREDP.ITYS. ID- IFDS. ISTAB. IPROP 
CCMMON/SPAR/MASSP.LW.XPP.YPP  »LB 
COMMCN/SURGE/SU( 7 ) 

COMMON/XSURC/XU . XUU . XUUU 
COMMON/SWAY/SW ( 15 ) 

COMMON/ YAW/ YA( U ) 

CCMMON/MTER/VCAR  .RHO.  ABS.CONI  ,CCN2 
CCMMON/RES I ST /VEL (AO). RES 1(40) 

COMMCN/VELE/UELC IOC) 

COMWON/POSTN/X1.Y1 .21 
COMMON/ GECM/AL . RW . G . AET . HW . HW 1 
CCMMSN/PRCP/ALE.P.EY.DI A.ANIU 
CCMmon.'CTnR/XC  <  99  ),YC(  99  ),ZC(  99  ).TC(  99) 

COMMON/ INT1/IC 
CCMMCN/DCC/UC .ALPHA 
CCMMON/UEPT/RLX .RLY -RL2 
COM“CN/SLCPE/PDRXX.PDRXY.PDRYX.PDRYY 
CCMMON/SLAN/RXX6 .RYY6 .RXX.RMY 

CPEN  (UNITs35.FILE*’ TANKER2 * .STATUS* ’OLD’ ) 

OPEN  (UNIT* 1. FILE* ‘RESIR1 .STATUS* 1  NEW ’  ) 

OPEN  (UNIT*2.FILE*,RES2R‘ .STATUS* ’NEW* ) 

OPEN  (UNIT*3.FILE*'RESJR' .ST ATUS* 1  NEW  ) 

OPEN  ( L'NIT*A  .F ILE*  'RESAR *  .STATUS* ‘NEW  ) 

CALL  INPUTDUO) 

VCAR  «VCAP«1 .689D0 

AL  *LW*LB«0.30A8D0 

ALE  *AL 

MAT2  *0 

IFLCWrl 

ILMAX* 1500 

EPS  *  1 .  D-5 


TOWOOOIO 
TCW00020 
TCW000S0 
TOWOOOAO 
TDW00050 
TCW00060 
TOWOOO  70 
TOW00080 
TOW00090 
TOWOOIOO 
TOWOOl 10 
TOW00120 
TOW00130 
TOWOOIAO 
TOWOOl 5  0 
TOW00160 
TOWOOl 70 
TOW00280 
TOWOOl 90 
TOW00200 
TOW00210 
TOW00220 
TOW00230 
TOW002A0 
TCW00250 
TCW00260 
TOW00270 
TCW00280 
TCW00290 
TOWOOJOO 
TOWOOJIO 
TOW00320 
TOW00330 
TOWOOJAO 
TOW00350 
TOW00360 
TCW00370 
TOWC0380 
TOW00390 
TOWOOAOO 
TOWOOAIO 
TOW00A20 
TOWOOAIO 
TCWOOAAO 
TOW00A50 
TOWOOA60 
TOW00A70 
T0W00A80 
TOW00A90 
TOW00500 
TOW00510 
TOW00520 
TOW00530 


WRITE 

(" 

1001 ) 

TOW005A0 

READ 

(» 

*  ) 

A1  ,A2 

TOW0C550 

WPITE 

(« 

1002  ) 

TOW00540 

READ 

(« 

"  ) 

NUM1 

TDW00570 

WRITE 

(• 

1003  ) 

TOM00580 

READ 

(• 

*) 

B1.B2 

TOW00590 

WRITE 

(« 

I00A) 

TOWC0600 

READ 

(« 

« ) 

NUM2 

TOWOOilO 

TOW0C420 

WRITE 

(• 

1005) 

TOW00430 

READ 

(« 

M 

1KB 

TQW006AJ) 
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WRITE  C**  - 1006)  TOW00650 

READ  (*•■)  IDS  TOWOO66O 

DO  l  1  =  1 .NUM1  TOW00670 

WRITE  (-.2001)  1 »NUM1  TOW0068Q 

LW  »A1*(A2-AI)*(I-I J/(NUM1-1)  TOM00690 

AL  »LW*LB*0 . 3 048 DO  TOW00700 

ALE*AL  TOW007I0 

XPP*B1  TOW00720 


CALL  ST AB I L  C I VV i VV » I SOL ) 

SET  V*VV(K)  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 
V*VV( 1KB) 

IF  CDABS(V).GT.l.DO)  STOP  1111 
CALL  ECUILBCV-X.RES.RX.RY) 

CALL  LINEAR (X.RE$.A.RX*RY) 

CALL  RGU  •  6  •  A .  WR  ,  WI  • MATZ*2 » I VI *SV2 * IER l ) 

IF  (IER1.NE.0)  STOP  2222 
IF  (IDS.EO.l)  CALL  DEGSTB C DECS • WR ) 

IF  (ICS.E0.2)  CALL  DS1CDE03.WR) 

DECSOO»DEO$ 

XPOO  «XPP 

L  «0 

DO  2  J*2.NUM2 

XPP»B1* (82-81  )«<J-1)/(NUM2-1) 

CALL  STASILdVV*  VV. I SOL ) 

SET  V*VV ( K )  FOR  BIFURCATION  ANALYSIS  OF  K-TH  EQUILIBRIUM 
V»VV( IKB  ) 

IF  CDABS(V).GT.l.DO)  STOP  Mil 
CALL  EQUILB(V.X.RES.RX.RY) 

CALL  LINEAR(X.REO.A.RX.RY) 

CALL  RG ( 6 • 6 • A . WR . WI . MATZ «Z« 1V1 • SV2 * IER I ) 

IF  (IERI.NE.0)  STOP  2222 
IF  (IDS.EO.l)  CALL  DEGSTB (CEOS. WR) 

IF  (IDS.E0.2)  CALL  DSl(DEOS.WR) 

DE0SNN.DE03 

XPNN*XPP 

PR«DECSOO*OEOSNN 

IF  (PR.GT.O.DO)  GO  TO  5 

L*LM 

IF  (L.GT.O  STOP  1000 
IL  =  0 

X-  0«XP00 

XPN*XPNN 

CECSO*DEOSOO 

DECSNrCEOSVN 

XPL»XPO 

XPR*XPN 

CEOSL*DEOSO 

DEOSR-DEOSN 

XPP*(XPL*XPR)/2,D0 

CALL  STABILC IVV»VV. I SOL ) 

V»VV( IKB) 

IF  (DABSCV).GT.l.DO)  STOP  1111 
CALL  EOUILB(V.X.REC.RX.RY) 

CALL  LINEAR(X.RES.A.RX.RY) 

CALL  RGUU.A.WR.WI.MATZ.Z.JVI *SV2.IER1) 

IF  (IER1.NE.0)  STOP  2222 
CALL  DEGSTBIDEOS.HR ) 

DEOSMtDEOS 
XPM*  XPP 


TOW00750 
TOW00740 
TOW00750 
TOW00760 
TOW00770 
TOW00780 
TOW00790 
TOW00800 
TOW00810 
TOW00820 
TOWO 08 20 
TOWO 0840 
TOW008S0 
TOWOO80O 
TOW008  7  0 
TOW 00880 
TOWO 08 90 
TOW00900 
TOW00910 
TOW00920 
TOW00930 
TOW00940 
TOWO 0950 
TCW00960 
TOW00970 
TOW00980 
TOW00990 
TOWOIOOO 
TOWOIOIO 
TCH01020 
TOWO l 030 
TOWO 1040 
TOWOIOSO 
TOW01060 
TOW01070 
TOW01090 
TOW0109C 
TOWOllOO 
TOWOUIO 
TOW0I120 
TOHOllSO 
TOWOIUO 
TOW01150 
TOWOIUO 
TOW01170 
TOWOIUO 
TOWOIUO 
TOW01200 
TOWO 1210 
TOW01220 
TOWO 1230 
TOWO 1240 
TOW012S0 
TOWO 1260 
TOWO 12 70 
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PRL*D£0GL"DE03M 

TOW01280 

PRR«DEC3R*CE0SM 

TOWO 1290 

IF  (PRL.GT . 0 . DO )  GO  TO  5 

TOWOI 500 

XPO=XPL 

TOW01510 

XPN=XPM 

TOWOI 520 

DE030*DE0SL 

TOWOI 550 

DE03N=DECSM 

TOWOI 540 

IL*IL*1 

TOWOI 550 

IF  (1L.GT.ILMAX)  STCP  5100 

TOWOI 560 

DIF»DABS(XPL-XPM) 

TOWOI 570 

IF  CD1F.GT.EPS)  GO  TO  6 

TOWOI 580 

XP«XPM 

TOWOI 590 

GO  TO  4 

TOW01400 

5 

IF  (PPR.GT.O.DO)  STCP  5200 

TOW0I410 

XPO«XPH 

TOWOI 420 

XPN=XPR 

TOW01420 

DE030«DE0SM 

TOW01440 

DEOSN=DEOSR 

TCW01450 

IL* 1L* 1 

TOW01460 

IF  (IL.GT.1LMAX)  STOP  5100 

TOWO 14  70 

DIF*DAB3(XPM-xPO ) 

TCW01480 

IF  CCIF.GT.EPS)  GO  TO  6 

TOW01490 

XP*XPM 

TCW01500 

4 

WRITE  (L-101  XP.LW 

TOW01510 

5 

XPOO=XPNN 

TOWOI 520 

DEOGOC  =  DEOSNN 

TOW01550 

2 

CONTINUE 

TOW01540 

1 

CONTINUE 

TOU01550 

STCP 

TOW01560 

10 

FORMAT  (2020.10) 

TOW01570 

1001 

FORMAT  (•  ENTER  RANGE  CF  Lw  VARIATION* ) 

TCW01580 

1002 

FORMAT  (*  ENTER  NUMBER  CF  INCREMENTS  IN  Lw*  ) 

TCW01590 

ioc: 

FORMAT  C  ENTER  RANGE  OF  Xp  VARIATION*) 

TOW01600 

1004 

FORMAT  (*  ENTER  NUMcSQ  OF  INCREMENTS  IN  Xp’) 

TCWC161C 

100S 

FORMAT  (*  ENTER  EC'JILIBRIUM  NUMBER*) 

TOW01620 

1006 

FORMAT  (*  ENTER  DEGREE  OF  STABILITY  CONTROL*) 

TOWC1650 

2CC1 

FORMAT  (215) 

TOW01640 

END 

TOWOI 650 

C 

TOW01660 

SUBROUTINE  DS1 (DEOS.NR ) 

TOW01670 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-2) 

TOWOI 680 

DIMENSION  NR(6) 

TCW01690 

DECS  1 *- 1 .d:o 

TOWOI 700 

DC  1  1*1.6 

TOWOI 710 

IF  (WR(I J.LT.DEOSl )  CO  TO  1 

TCW01720 

CEOS  1  *wR { I ) 

TOWOI 750 

1J*I 

TOWOI 740 

1 

CONTINUE 

TOWOI 750 

DEOS2*-l .D20 

TOWOI 760 

DO  2  1*1.6 

TOWOI 770 

IF  (IJ.EO.I)  GO  TC  2 

TOWOI 780 

IF  (WR(  I  KLT.DECS:  )  GO  TO  2 

TOWOI 790 

DEOS2*WR( I ) 

TOW01800 

IJJ»I 

TOW0181 0 

2 

CONT  I  N’JE 

TOW01820 

DEOS*-l .050 

TOW 01850 

DO  5  1*1.6 

TOW01840 

IF  (I.E0.IJ.0R.1.E3.UJ)  GO  TO  5 

TOW01850 

IF  ( WR  ( I ) . G£ . DECS )  DECS* WR ( I ) 

TOWO I860 

5 

CONTINUE 

TOU01870 

RETURN 

TOWO 1880 

END 

TOW01890 
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